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NUMERICAL MODELING OP STIMULATION OP INDUCED 
MAGNETOSPHERE DURING INTERACTION OP SOLAR WIND 
WITH THE lONOSPI-IERE OP VENUS 

A. S. Lipatov 


INTRODUCTION 

Measurements of plasma and magnetic field in the vicinity 
of Mars [1-4] and Venus [4-6] have discerned a "head” non- 
colliding shock wave, a weak magnetic force on Mars - 50 
and the absence of a field on the surface of Venus. 

At the present time, gasdynamic models of atmospheric 
flov7 of Mars and Venus [7-8] are known. It has been supposed 
that the external boundary of the atmosphere is a tangential 
break, and the braking pressure in the transitional sphere is 
balanced by the pressure of the atmosphere. 

Also known are quasi-hydrodynamlc evaluations of the value 
of induced magnetic field associated 'with the diamagnetism of 
photoions and conductivity of the ionosphere, which do not con- 
sider reverse influence of the field on flew in the transitional 
sphere [9]. 

Direct numerical modeling of the Interaction of solar wind 
plasma with the ionosphere [10 J has shown that in the case of 
subsonic flow, an induced magnetosphere may exist. Hovj'ever, 
further calculations on supersonic regimes have shown that 
plasma and solar wind currents are locked on the_ ionosphere. 


^Numbers in margin Indicate pagination of original foreign text. 
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This fact has led to the necessity of considering inter- 
action processes in the ionopause in the overall flow diagram. 

This work is devoted to the study of: 

1) electrodynamic processes in the ionopause; 

2) the structure of the induced magnetosphere. 




MATHEMATICAL FORMULATION- OP THE PROBLEM 


/I 


Solar wind model . Proceeding from solar wind conditions 

we will describe electrons by equations of magnetic hydrodynam' 
ics, and ions by function of distribution in a dimensionless 
form; 


generalized Ohm's law 
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equation for ions 
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Maxwell's induction equation 

where V = < \4 > — average macroscopic speed o,f’ plasms, 

A/^ f A/i — al venous and sonic Mach numbers, W cp t.UJi 
— unstimulated speed of solar wind and ion gyro frequency ,• 

R* — characteristic dimension of the problem — outside radius 

of the ionosphere, _ 4'^ V/co Co — Reynold’s magnetic 

^ ■ C ® 

number, — effective adiabatic indicator (in computations 
^ = 5/3), — conductivity. 


The transition from dimensional variables to dimension- 
less is determined by formulas: 

■ ^ i' ■ X- n* ; 
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E - B‘ ; P= P^P' 


Ionosphere model . Let us suppose that the parameters of 
the planet’s atmosphere are known. Then the dynamics of the 
ionosphere is described by equations of .motion, field, and in- 
dissolubility: ' 
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— concentration of charged and neutral particles, 

— coefficients of photoionization and recombination, 

— speed and pressure in the ionosphere, 

— force of friction with neutral particles. 


The presented mathematical model of the ionosphere is 
rather complex. Since in this work we are primarily Interested 
in processes in the transitional sphere and- in the ionopause, 
we will examine a simplified model of the ionosphere: 


1) The distribution of charged particles in the iono- 
sphere is stationary; 

2) convection in the ionosphere in the first approxi- 
mation is absent ^ = 0 ; 

3) "throw-off" of magnetic field from the ionosphere 
is determined by effective conductivity 

lol R.en 

— Reynolds number determined by Coulomb collisions); 

^1) the ionosphere has the form of a spherical ring with 
thickness • ’ 


Initial and boundary conditions . It is assumed that at an 
initial moment in time, the ionosphere is submerged in a homoge- 
neous interplanetary magnetic field. The parameters of plasma 


in the vicinity of the ionosphere are also unstimulated. At 
the outside boundary of the computed sphere (half the spheri- 
cal ring, Pig. 1) = Ha , parameters of the oncoming 

plasma flow and fields correspond to interplanetary, and 
Wo* i boundary sphere Q ^ /2 > homogeneous 

boundary conditions of type II 2. /-/ - are set for the mag- 

*<) 0 

netlc field. At the Inside boundary of the ionosphere 

= I - the magnetic field is considered unstimulated. 

II. NUMERICAL MODELING METHOD 


Discrete model of plasma . The idea of the "particle in a 
cell" method belongs to Harlow and has been used repeatedly 
for modeling not only laboratory (see, for example, [ll-l4j ) , 
but also cosmic plasma [10, 15-17]. This work utilizes the 
"CIC" ("clouds in a cell") method for computing the ion com- 
ponent of solar wind. The oncoming flow is approximated by a 
cold bundle of ions having spatial structure. Since equations 
of the field (1,2,6, 7) are solved in planes ^ = P, /2 
(solution for 0 < is Interpolated by trigonometric 

series), it has turned out to be expedient to inject plasma ac 
cording to the following diagram. Pigs. 2, 3. This Initial 
distribution provides full approximation of parameters of 
plasma in sections y s 0^ ^ /2 with a minimal number of 
macroparticles /V^= 100 X 15. X 71 

Differential method of solving Maxwell’s' equations . In 
solving field equations, decomposition into Pourler series is 
utilized by i 

^ ^ CO.!. 2 9 
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* ^ 2 . — value of field with ~ 0 and 
The induction equation in planes _y* = 0 and ^ = JL /2 
is as follows ; 


hi = L Lki.-I-^H'^J ' y^/ "’ ^ 


Li Li = v/k A, ^ Li ^•) . 

4 L-i = Tkyy > L-i ' 

In solving (9), a uniform grid according to R and ^ (.30 x 

20), Pig. 4, is used. The Krapka-Nicholson method is used 
cd = 0,5) , differential operators /j ^ and ^ have the approxi- 
mation^ ( ^nd Q ) > respectively. Mixed derivatives 

are approximated by evident differentials. 

Limitational conditions ^ l\9= 0 kst • ks> • 

are placed at axis ^ for k'. 


In calculating the dynamics of ions and their concentration 
grid II is Introduced, shifted by A -i /2. .-it /2, ■^6>/2 


relative to Grid I, Pig. i|. 


The equation of motion fs approximated by a non-evident 
method with approximation Q ), 

III. RESULTS OP CALCULATIONS 

Calculations are performed for the following parameters 
of oncoming flow: 

M^=M3 =10; = 0,086; 36 =5/3; 

l2u = 10 ^ - 10 ^; 
r I»5o. 


As a result of oncoming flow of plasma onto the ionosphere, 
electrons and ions penetrate to the depth of gyro-radius order 
'tc e j since the electrical fiel'd in the ionosphere is small 

Eu Hic . 

Since with initial magnetic field, the dimensionless Larmor 
radius of ions <’ = Wp . l — = 0.0 86 Is compared with the 

thickness of the ionosphere du =0.1, ions, forming a current 
layer, fill up the ionosphere. Electrons having 'Ye. 3 grades 
lower than ^ practically do not penetrate into the ionosphere. 
The final Larmor ion radius causes viscous interaction of plasma 
in the transitional sphere with the ionosphere. This fact 
significantly distinguishes this interaction model from gas- 
dynamic models 17-8]. 

Let us consider in detail the flow picture in the regime 

» yiu. - 10^ moment in time T-* 0.85. Pigures 5--6 ‘ 

give tne distribution of solar wind plasma concentration in the 


transitional sphere and the ionosphere in an equatorial and 
meridional plane. Evident from these graphs is a sharp in- 
crease in concentration directly behind the front of the 
shock wave. Away from the axis with , stimulation 

of concentration drops. 

Figures 7-8 give the distribution of magnetic field 
module in an equatorial plane and the configuration of mag- 
netic force lines in a meridianal plane. Directly behind 
the front of the shock wave, the magnetic field Increases ■ 
several times over, following the concentration. A conden- 
sation of force lines is observed near the ionopause, the 
plasma "frozen in" is disrupted and the magnetic field does 
not correspond to concentration. Pinal conductivity causes 
diffusions of the magnetic field into the ionosphere. The 
characteristic thickness of the "stagnant" sphere with strong 
magnetic field may be compared with the gyroradius in the 
transitional sphere. 

Figures 9-10 give plasma current lines in the equatorial 
and meridianal planes. It is evident from the graphs that the 
current lines bend around the zone of the strong magnetic field. 

In spheres of relatively weak magnetic field, plasma flows into 
the ionopause, partially Including the interplanetary field. 

Prom calculation results, we can see the asymmetry in the 
flow picture in the equatorial and meridianal planes. Moreover, 
due to the Influence of the final Larmor ion radius, a displace- 
ment of zero current line in the equatorial plane was observed. 

Figure 11 gives the dependence of maximal strengthening /lO 

of the magnetic field in the ionopause from Reynolds magnetic 
number in the ionosphere . It is evident that with re- 
duction of , the diffusion of the magnetic field through 
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the ionosphere grows, and the ''piimplng" of the field in the 
ionopause is reduced. It should be noted that speed of 
convection in the ionosphere corresponding to effective diffu- 
sion may be evaluated as ]/u 'v ^ 

— characteristic diffusion scale). 

Figure 12 gives the change in thickness of the transi- 
tional sphere in the subsolar point in time depending on 

. Calculations have shown delay in formation of re- 
flected particles and shock wave due to strong diffusion of 
the field in the ionosphere, and consequently a large Larmor 
ion radius. The thickness of the transitional sphere in our 
model exceed?^ the value of gasdynamic models [7-8] . This re- 
sult is associated primarily with the rejection of the assump- 
tion that the ionopause is a tangential break, and consequently 
with >increase of viscosity effects. 


CONCLUSIONS /II 

1. During flow around the ionosphere, a considerable mag- 
netic field is stimulated in the i’onopause, associated with 
the final size of Larmor ion radius. 

2. The amplitude of the magnetic field in the ionopause 
is significantly determined by the mechanism of "transfer" of 
the field through the ionosphere. Particularly, an accounting 
of real convection of the photoionosphere i^ necessary for 
evaluation of parameters of plasma and field in the ionopause. 

3. The magnetic hydrodynamic boundary layer stimulated . , 

in the ionopause determines the flow in the transitional sphere, 

4. Nonhydrodynamic approximation has additional asymmetry 

I 

of flow relative to the meridional plane. 
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FIGURE CAPTIONS 


Fig. 1. 
Fig. 2. 
Fig. 3. 
Fig. 

Pig. 5. 
Pig. 6. 
Pig. 7. 
Pig. 8. 
Fig. 9. 
Pig. 10. 
Pig. 11. 

Pig. 12. 


Geojnetry of computational sphere of the problem 
Distribution of particles in planes X - Y 
Distribution of particles in planes X-Z, Y-Z 
Differential grids for computing field -k — and 
plasma parameters — . — 

Distribution of concentration in plane X-Z 
Distribution of concentration in plane Y-Z 
Dines of levels H in plane X-Z 
Magnetic field force lines in plane Y-Z 
Current lines in plane X-Z 
Current lines in plane Y-Z 

Dependence of value of magnetic field in the iono- 
pause on Reynolds magnetic number ■' the ionosphere 
Change in thickness of the transitional sphere 
through time for various Reynolds numbers of ..the. . 
ionosphere 
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